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Abstract 

Q! We study infrared conformality of the twelve-flavor QCD on the lattice. Utilizing the highly 

(D. 

■ improved staggered quarks (HISQ) type action which is useful to study the continuum physics, 

we analyze the lattice data of the mass and the decay constant of a pseudoscalar meson and the 
mass of a vector meson as well at several values of lattice spacing and fermion mass. Our result is 
consistent with the conformal hypothesis for the mass anomalous dimension 7^ ^ 0.4 — 0.5. 
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I. INTRODUCTION 



There has been a renewed interest in the study of the QCD with a large number of the 
massless fermions in fundamental representation ("large Nf QCD") in the context of walking 
technicolor having approximate scale invariance and large anomalous dimension 7^ ~ 1 l| 
The model was proposed to cure the fatal flavor- changing neutral current (FCNC) problem of 
the original technicolor and predicted a technidilaton as a pseudo-Nambu-Goldstone (NG) 
boson of the approximate scale invariancci, based on the ladder Schwinger-Dyson (SD) 
equation with nonrunning (scale-invariant) gauge coupling. (See, for a review, Ref. [t].) 
Such an approximate scale-invariant dynamics may in fact be realized in the large Nf QCD: 
The perturbative two-loop beta function predicts a nontrivial infrared (IR) fixed point 
(0 < a* < 00) in the range of 9 < A^/ < 16 in the asymptotically free SU(3) gauge 
theory j^, 9], where Nf is the number of massless flavor. The theory has an intrinsic scale 
A, analogue of Aqcd, which characterizes scale-symmetry breaking (scale anomaly) in the 
ultraviolet (UV) region (/i > A) where the coupling runs as an asymptotically free theory, 
a{fi) ~ l/log(yu/A), in the same way as the ordinary QCD. However, the theory in the IR 
region (0 < // < A) governed by the IR fixed point respects approximate scale invariance 
due to the almost nonrunning (walking) coupling (IR conformality) , a{fi) ~ a^,. If such 
an IR fixed point survives the nonperturbative effects, it would imply that the theory is in 
the deconfined and chiral-symmetric phase, the so-called "conformal window" (though the 
theory has no conformality in the UV region), which is expected to exist for a certain range 
of large A'^^. 

Actually, the conformal window and phase structure of the large Nf QCD was studied by 
the (improved) ladder SD equation in a way to simply replace the coupling by the two-loop 



running coupling having the IR fixed point mentioned above |lOl . Ill| : For 9 < Nf < 11, 
the IR fixed-point coupling exceeds the critical coupling a* > acr and hence triggers the 
spontaneous chiral symmetry breaking (S^SB) to give rise to the dynamical generation of 
the fermion mass rrioi^ 0), which actually washes out the would-be IR fixed point, since the 
fermions, once acquiring mass, are decoupled from the beta function in the very IR region 
below the scale of m^- The conformal window thus is expected to be 12 < Nf < 16. If it 
^ For subsequent similar works without notion of scale invariance and anomalous dimension see 2, 

4|. 



31. For 



earlier work on this line based on purely numerical analysis see 
^ The technidilaton may be identified with 125-126 GeV boson observed recently at LHC 

discussions can be seen in Ref. [6[ which appeared after submission of this paper. 



5| . More detailed 



is the case, the walking technicolor should be realized in a slightly broken phase very close 
to the conformal window, Nf ^ 12, such that nir) <^ h., with m/) being on the electroweak 
scale order and A being usually identified with the extended technicolor scale. Although the 
IR fixed point actually disappears at ~ m£)(^ A), the coupling is still almost nonrunning 
as a remnant of the would-be IR fixed point for the wide IR region rriD < fi < A, which is 
relevant to the physics with the critical value 7^ — 1 valid all the way up to the scale A. 

In the conformal phase with the bare fermion mass rrif, the low-energy behavior of the 
hadron spectra should obey the same scaling relation (hyperscaling relation) near the IR 
fixed point, M ~ my^^'^'^'*\ where M is the hadron mass and 7* is the mass anomalous 
dimension 7^ at the IR fixed point. For small rrif, we can expect that 7^ — 7* in the 
wide region all the way up to the scale A where the IR conformality is operative. On the 
other hand, if the chiral symmetry is spontaneously broken (not in the conformal phase), 
the low-energy physics could be described by the chiral perturbation theory (ChPT). 

As a powerful nonperturbative study of such a problem, one can use lattice QCD sim- 
ulations, which can, in principle, determine the phase structure of the SU(3) gauge theory 
with a various number of fermions through investigating nonperturbative running of gauge 
coupling, low-energy spectra, chiral condensates, and so on. Furthermore, when the theory 
is inside the conformal window, one can also calculate critical exponents such as the mass 
anomalous dimension. A similar analysis could be approximately applied to the slightly 
broken phase (m^) ^ ^ A) refiecting the remnant of the conformality, which is vital to 
the study of the walking technicolor. 

In lattice calculation, there exist effects of the finite size boundary and finite fermion 
masses, which breaks the IR conformality explicitly and also deforms the ChPT. In the case 
of the conformal phase, it may cause serious uncertainties, since the IR scale is given by 
the inverse lattice size or input nif but not by the intrinsic scale A(^ rrif^L"^) which 
behaves as the UV scale where the IR conformality breaks down, in sharp contrast to the 
ordinary QCD where the IR scale is given by the intrinsic scale Aqcd{^ i^f, L~^). Actually, 
the finite mass mj as well as finite size L can easily distort the hyperscaling relation, while 
finite mj can also trigger the dynamical mass m/) of the fermion since the coupling should 
blow up due to the decoupling of the fermion in the IR region fiKnij, thus mocking up the 
ChPT. 

On the other hand, even in the chiral-broken phase, there still can exist an approximate 
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hyperscaling relation, as far as rriD ^ rraj ^ A, as a remnant of the approximate scale 
invariance, while there exist finite size effects and higher-order contributions which make 
it difficult to understand the ChPT particularly in the large A'^^ QCD. Here, we should 
point out that the validity region of the ChPT is extremely restricted near the conformal 
window where the decay constant of vr (/tt) in the chiral limit is expected to be very small 
and, particularly, Nf is large. For the consistency of the ChPT 12|, it is required that the 

is]- 15 1 for the mass of vr (M^) is small, 

iV/f^V<i (1) 



expansion parameter 



which, however, could become easily violated when simulation is made for away from 
chiral limit. 

There are many lattice results on the large Nf QCD in recent years. In addition to the 



pioneered works 



16 



181 ] for the study of the phase structures, there has been a lot of progress 



on the lattice calculations such as the renormalized running coupling, hadron spectra, finite 



temperature transitions, etc. 



19|-[30|]. In particular, for the A'^^ = 12 the system has been 



widely investigated by the lattice approach focusing on determining the phase structure: A 



recent 
phase 



arge- volume analysis for the hadron spectra gives the result favoring the chiral-broken 



24l |. while other groups with similar analyses using the same data conclude that this 



model favors the conformal phase 3l|, . Thus, it is controversial at present whether this 



theory is in the conformal or the chiral-broken phase. 

Another concern is the value of /3(= Q/g"^) in simulations. It is suggested that there 
are some phase structures inherent to the lattice model even for a small number of flavors. 
Actually, following earlier suggestions Jq, 1i| [. there have been some studies about the lattice 



phase diagram in the large Nj QCD 



21 



29 



33| which, in fact, indicate a nontrivial lattice 



phase with a bulk transition at /3 other than the critical value of /3 which is associated with 
the IR fixed point we are interested in. Thus, it is important to survey the (3 dependence 
for the investigation of the low-energy physics. 

In this work, we study the phase structure of the twelve-flavor QCD. Preliminary results 
were given in Ref. 3^. We utilize the highly improved staggered quarks (HISQ) type ac- 



tion 



35 



,136| to reduce the discretization error. This is the flrst attempt to use the HISQ-type 
action for the large A'^^ QCD. A salient feature of our work is that we perform simulations of 



N 



f 



12 in comparison with other cases of Nf = (0), 4, 8, 16 which we do simultaneously with 



the same systematics based on the same HISQ-type action SJ, l37|. We make simulations 
at two values of /3 (= Q/ g"^) (or two lattice-spacings), with g being the bare gauge coupling, 
to study the lattice spacing dependence of the data for the reason we described above. Our 
results suggest that two coupling regions are consistent to be in the asymptotically free 
region, where we discuss lattice spacing dependence of the physical quantities. 

We investigate several bound-state masses such as the pseudoscalar meson vr (correspond- 
ing to the exact NG boson if it is in the broken phase) and vector meson p as well as the 
decay constant of vr, by varying the fermion bare mass mj. 

Concerning the controversy about the previous lattice studies of Nf = 12 QCD mentioned 
before, it should be pointed out that there is a problem in judging whether this theory is 
in the conformal or the chiral-broken phase from the fit analysis using lattice data of the 
hadron spectra with finite fermion bare mass rrif, which fully depends on the fitting forms 
of the ChPT fit in the chiral-broken phase as well as on the hyperscaling relation in the 
conformal hypothesis. Indeed, we do not know the systematics of the corrections due to 
the finiteness of the mass trif and the volume size L in either case, which could make some 
confusion. 

In the present work, we introduce an alternative quantity for the scaling test of the 
conformal hypothesis with finite volume in the analyses of the simulations for Nf = 12 QCD. 
Using this method, it is possible to analyze the data without any specific function form as a 
fitting function to the data. We can discuss possible finite mass and finite size corrections. 
We find that the results are consistent with hyperscaling for 7* ~ 0.4 — 0.5. This can be 
compared with the value 7^^ ^ 0.80, the value calculated through the ladder Schwinger- 
Dyson equation 7^° ~ 1 — ^1 — ^ [l|, [ssl for two- loop a^, ^ 0.754 and acr = 7r/4 [fj. 
Moreover, for concrete understanding of possible corrections, we try to perform several fits 
with some correction terms. 

We also perform the ChPT analysis of our data. It turns out that our data are far from 
satisfying the consistency condition in Eq. ([T]). 

Our paper is organized as follows. In Sec. [Tll we explain our model and simulation setup. 
The results of the numerical simulations for the spectra are shown in Sec. IIIII Further study 
of dimensionless ratios composed of these measurements will be performed for the study 
^ Using three-loop (four-loop) a* in the MS-bar scheme in this expression, one obtains 7^^ ~ 0.33 (0.37). 



Note that 7^^-* can be written as 7^^-* ~ 1 — yl — 27^"'°°^. This is also compared with perturbative 

two-loop (7* ~ 0.77), three-loop (7* ~ 0.31), and four-loop (7* ~ 0.25) calculations 
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of the hyperscaling. The consistency of the asymptotically free domain is also discussed 
from the /3 dependence of these observations. Our main result is given in Sec. llVt where 
we introduce a method for the scaling test for conformal hypothesis and analyze the data 
using this method. We deduce the anomalous dimension from our analysis. In Sec. |V] we 
perform the analysis based on the fit function which includes the corrections to hyperscaling 
as suggested by certain models. In Sec. IVIt we briefly summarize the results on the analysis 
of the SxSB scenario. The summary and discussion are given in Sec. I VI II In Appendix |Xl 
we show the results for other spectra than those discussed in the text. In Appendix [Bl we 
show the numerical details of the fit results on the finite size hyperscaling. In Appendix O, 
the details of the fit analysis based on the ChPT are shown. 

II. THE BASIS METHODS FOR THE STUDY OF THE CHIRAL PROPERTIES 
A. Continuum theory 

Our target is the SU(3) gauge theory with Nj = 12 of massless Dirac fermions in the 
fundamental representation. Investigating the spectral quantities in mesonic channel of 
the mass-deformed theory, we try to determine the phase and the quantities which are 
characteristic to the associated phase. We use the technique of lattice gauge theory, whose 
continuum counterpart is briefly described here. The continuum Euclidean action reads. 



where ifjq denotes the 12-flavor fundamental fermions with degenerate mass ruf. The flavor 
nonsinglet bilinear operators = tp'-f^a'^tp for the pseudoscalar and = ■ip'-f4'-fia"'ip for the 
vector channel are used to interpolate the bound state with the particular quantum numbers 
under flavor and Lorentz symmetries. Here, a" is a generator of the SU(12) flavor symmetry 
group. The ground-state mass, Mh where /J = vr for the pseudoscalar (O = P") or H = p 
for vector (O = V^) characterizes the asymptotic fall off Go ~ e~*^^* of the Euclidean 
correlation function with the zero-momentum projection 






(3) 
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calculated using the action Eq. (|2]). The pseudoscalar decay constant is obtained through 
the matrix element of the pseudoscalar operator, 



TD £ — » 

F. = jf,{o\P^{o)K;p = o), (4) 

using partially conserved axial current (PCAC) relatione. 

If the theory is in the conformal window, Mh and obey the conformal hyperscaling 

Mh oc mf^, oc mf^, (5) 

where 7* denotes the mass anomalous dimension at the IR fixed point. This relation is the 
fundamental scaling appearing in the critical phenomena of the statistical mechanics. In the 



context of the conformal window in the large Nf QCD, see, for example, Ref. 40|. On the 



other hand, if the theory is in the phase of S^SB, leading mass dependence will be 

oc m/, F^ = co + cinif, (6) 

with Co 7^ 0, and the vector meson mass does not vanish in the chiral limit. 

The spectra obtained in our lattice simulation will be tested against these two hypothesis 
in the following sections. 



B. Lattice setup 

A staggered fermion formulation is used to define the lattice version of the Nf = 12 
SU(3) gauge theory. A theory with three degenerate staggered fields with the mass rrif has 
the degree of freedom of the A'^^ = 12 Dirac fermions and matches with the theories with 
twelve degenerate flavors in the continuum limit from the asymptotic free regime. Note 
that there is no problem of locality in this lattice theory, as the rooting trick to tune the 
number of flavors to Nf 7^ 4n theory is not used here. At a nonzero lattice spacing, where 
one can perform numerical computations, the 12-fold degeneracy does not hold, in general. 
This nondegeneracy would manifest itself as the nondegeneracy of the would-be degenerate 
12 X 12 — 1 mesons in the continuum theory. As the number of the effective light degrees 
of freedom is important to the critical phenomena, such as the IR conformality near the IR 
fixed point, this nondegeneracy needs to be made much smaller than other characteristic 

^ We use the convention as = v^/tt, where /^r = 93[MeV] in the real-hfe QCD. 
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scales of the theory. This can be achieved by taking the continuum hmit or using the lattice 
action which suppresses these effects at a given lattice spacing. 



We use a version of the HISQ 35|, |36|] to simulate this system. This action has the best 



continuum scaling among the staggered actions used so far in the real QCD simulations. 
For practical reasons, the tree-level improved Symanzik gauge action instead of the 1-loop 
tadpole improved one is adopted. Exactly the same setting, called HISQ/tree, but with the 
rooting for 2+1 flavor simulation, has been used for the QCD thermodynamics and has lead 
to an expectedly good reduction of the flavor-symmetry violation in the pseudoscalar mass 



splitting 4l|. We actually observe excellent results for the flavor symmetry in our Nj = 12 
simulations, which is briefly shown in Appendix lAl 

We use three lattice volumes (L,T) = (18, 24), (24, 32), and (30,40) and two lattice 
spacings with /3 = 3.7 and 4.0, where L and T are spatial and temporal length of the 
flnite lattice and /3 = Q/g"^ with bare gauge coupling g. Note that the aspect ratio T/L 
is kept flxed = 4/3 while the volume is changed. In this way, we only have one IR scale 
1/L, which is ideal for the flnite-size scaling analysis. For the quark mass rrif, we take 
various values: mj = 0.04,0.05,0.080,0.1,0.12,0.16,0.2,0.24, where 0.04 and 0.24 are only 
for /3 = 3.7 and 4.0, respectively. The gauge conflgurations are generated by the hybrid 
Monte Carlo algorithm using MILC code ver. 7 with some modiflcations to suit our needs. 
The boundary condition is set periodic for all except temporal boundary for fermions, which 
is made antiperiodic. We accumulate 400-1200 thermalized trajectories for each ensemble. 
Calculation of the mesonic correlation functions is performed at every 2-10 trajectories, thus, 
we have 70-400 samples for each ensemble. The error analysis is performed with the standard 
jackknife method with a suitable bin size 10 — 100 depending on the ensemble parameter. 

For comparison, we also analyze the Nf = 4 theory which is in the phase of S^SB, where 



the same lattice action is used. Some detail of the Nf = 4 simulation is given in Ref. 37 1 
and the forthcoming publication. 



III. NUMERICAL RESULTS AND PRIMARY ANALYSIS 

We measure the two-point correlation functions of the staggered bilinear pseudoscalar 
operator which corresponds to the NG mode associated with the exact chiral symmetry of 
the staggered fermions. The corresponding spin-flavor structure is (75 ® .^5), denoted by 
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"PS" in Ref. 42|]. The random wall source is used for the quark operator for the bilinear, 



which becomes a noisy estimator of the point bilinear operator with spatial sum at a given 



time slice to- We combine quark propagators solved wit 



1 periodic and antiperiodic boundary 



conditions in the temporal direction (see, e.g., Ref. 43|). In this well-known technique, the 
temporal size is effectively doubled, which enables us to have sufficient range for the fitting, 
which only takes into account the ground state. Denoting such a vr correlator as Cps(t), 
we use C'ps(2t) = Cps(2t)/2 + Cps(2t - l)/4 + Cps(2t + l)/4 for the analysis. This linear 
combination kills the constant oscillation mode, which could originate from the single quark 
line wrapping around the antiperiodic temporal boundary. The results of effective mass 
calculated with C*ps(2t) for L = 30 are shown in Fig. [T] (see those getting plateau from 
above). Effective mass of the same vr correlation function with the Coulomb-gauge-fixed 
wall source is also plotted (ones getting plateau from below). Note that the plateau starts 
early enough to isolate the ground state even if the temporal size had been T = 24 (our 
smallest). The vr mass is obtained by the fit of the two-point correlators of Cps from a 
random source with double period by a fit function 

C'ps(2t) = C(e-^^-2i ^ g-Af.(2T-2t))^ 

where is a mass of the vr and the fit range is taken [tmini T]. Similarly we can obtain 
from this operator. The values of tmin are taken as tmin = 16 — 22 and 18 — 22 at /9 = 3.7 
and 4, respectively. 

We rneasure Mp from the staggered vector meson operator (7^74 ® ^i^i), denoted by PV 



in Ref. 42|]. The asymptotic form of the PV correlator may be written as 

G'pv(t) = C^{e-''''' + e-^^''(2T-t)) _^ c'2(-l)*(e-^^'"i* + e-^^^i^^r-t)) 

where M^^ corresponds to the mass of the axial vector meson. Since there exists a constant 
mode due to the wrapping- around effect, we use Gpy{2t) = Gpy{2t)/2 + G'pv(2t — l)/4 -(- 
Gpv(2t + l)/4. Again, the results of effective mass calculated with Gpv for L = 30 are 
shown in Fig. [2l The fitting range for the PV correlator is [tmin, T] with tmin = 10 — 12, and 
12 — 14 at /3 = 3.7 and 4, respectively, which is enough to isolate the ground state, i.e. p. 
We obtain Mp by the fit of the two-point correlators of Gpv by a simple fit function: 

Gpv(2t) = C{e-''^'' + e-*^^p(2^-2t))_ (g) 



All the raw results of Nf = 12 theory which are used in the next sections are shown in 
the Tables Hl lVIl In these tables, we also show the number of the thermalized trajectory for 
each parameter, by which we measure the masses and the decay constant. 

Besides these main channels, we study the masses of mesons made with local operators, 
a non-NG channel (7574 ® ^5^4) denoted by "SC" , and a vector meson (7^ ® ^j) denoted by 



"VT" in Ref. 



42|. They are compared against PS and PV in Appendix |Al which indicates 



the flavor-breaking effect of the HISQ is very small. 
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FIG. 1. Effective mass of vr at /3 = 3.7 (Left) and /3 = 4 (Right) on {L,T) = (30,40) with Cps 
using two types of quark sources. See the text. 
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FIG. 2. Effective mass of the PV correlator (Gpv) at /3 = 3.7 (Left) and /3 = 4 (Right) on 
(L,T) = (30,40). 



For a primary analysis, dimensionless ratios composed of these measurements will be 
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TABLE I. The results of the spectra on V = 
18^ X 24 at /3 = 3.7. iVtrj means the total num- 
ber of thermalized trajectories. 



TABLE n. The results of the spectra on V 
18^ X 24 at /3 = 4. 









Mp 




0.04 


1000 


0.3387(62) 


0.4614(85) 


0.0633(24) 


0.05 


1000 


0.3830(38) 


0.4899(71) 


0.0750(22) 


0.06 


1200 


0.4100(51) 


0.5144(49) 


0.0814(12) 


0.08 


1200 


0.4922(17) 


0.6139(33) 


0.1027(11) 


0.1 


700 


0.5717(17) 


0.7084(33) 


0.1185(6) 


0.12 


700 


0.6485(16) 


0.7905(38) 


0.1333(11) 


0.16 


700 


0.7878(13) 


0.9546(20) 


0.1590(9) 


0.2 


700 


0.9209(10) 


1.1041(32) 


0.1838(8) 



F^ 



0.05 1000 0.3993(67 

0.06 1200 0.4101(41 

0.08 1200 0.4551(22 

0.1 1200 0.5219(43 

0.12 1200 0.5960(18 

0.16 1000 0.7309(21 

0.2 700 0.8567(11 

0.24 700 0.9760(17 



0.4529(88) 
0.4687(115) 
0.5438(44) 
0.6133(89) 
0.6955(40) 
0.8564(62) 
0.9903(38) 
1.1225(29) 



0.0573(9) 
0.0674(11) 
0.0878(10) 
0.1009(19) 
0.1149(10) 
0.1378(11) 
0.1586(12) 
0.1785(14) 



TABLE III. The results of the spectra on V = TABLE IV. The results of the spectra on V = 



24^ X 32 at ;5 = 3.7. 






24^ X 


32 at /3 = 4. 






ruf 






Mp 


F^ 


ruf 






M, 


F^ 


0.04 


600 


0.3054(31) 


0.3569(86) 


0.0621(15) 


0.05 


1000 


0.3290(25) 


0.3819(81) 


0.0615(7) 


0.05 


600 


0.3549(16) 


0.4377(44) 


0.0750(7) 


0.06 


1000 


0.3677(13) 


0.4353(43) 


0.0714(7) 


0.06 


1000 


0.3986(22) 


0.4932(54) 


0.0835(7) 


0.08 


1000 


0.4459(11) 


0.5257(31) 


0.0875(9) 


0.08 


800 


0.4880(8) 


0.5991(25) 


0.1014(4) 


0.1 


1000 


0.5210(7) 


0.6165(27) 


0.1014(4) 


0.1 


800 


0.5683(17) 


0.6971(32) 


0.1175(6) 


0.12 


1000 


0.5946(11) 


0.6999(21) 


0.1147(5) 


0.12 


700 


0.6459(5) 


0.7891(24) 


0.1325(4) 


0.16 


700 


0.7308(9) 


0.8519(27) 


0.1378(5) 


0.16 


700 


0.7879(6) 


0.9529(22) 


0.1591(8) 


0.2 


700 


0.8557(5) 


0.9893(29) 


0.1579(4) 


0.2 


700 


0.9193(6) 


1.1031(21) 


0.1818(8) 













plotted against the tt mass. In the rest of this section, we write lattice spacing a explicitly 
and consider M^r as dimensionful quantity. As a reference, we have calculated the Nf = 4 
case as plotted in Fig. [3l which resembles a familiar QCD-type behavior. The upward 
tendency of F^/M^ toward the chiral limit is consistent with the S^SB [Eq. ([6])]. 
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TABLE V. The results of the spectra on V = TABLE VL The results of the spectra on 
30^ X 40 at /3 = 3.7. F = 30^ x 40 at /3 = 4. 



ruf 


AT 


1\ /F 


7\ r 
Mp 


77' 


ITlf 








0.04 


onr\ 

oOU 


0.3Uzo(14j 


0.3713(39) 


0.0D37(Dj 


0.05 


600 


0.3167(27) 


0.3671(56) 0.0634(8) 


0.05 


700 


0.3504(11) 


0.4302(20) 


0.0741(5) 


0.06 


700 0.3648(14) 0.4357(17) 0.0732(6) 


0.06 


600 


0.3990(15) 


0.4864(28) 


0.0835(5) 


0.08 


600 


0.4499(8) 


0.5301(13) 0.0901(5) 


0.08 


500 


0.4869(8) 


0.5949(15) 


0.1017(5) 


0.1 


600 


0.5243(7) 


0.6150(16) 0.1027(6) 


0.1 


500 


0.5670(7) 
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FIG. 3. Dhnensionless ratio F-j^/Mt^ as a function of aAf^r for Nf = 4 at /3 = 3.7. Due to the SxSB, 
the ratio diverges in the chiral limit. 

The left panel of Fig. H] plots the same quantity for Nf = 12 at the largest two volumes 
and two bare gauge couplings /3 = 3.7 and 4. If we look at the results for /3 = 3.7 (filled 
symbol), Ft^/M^^ tends to be flat for smaller vr masses, which shows clear contrast to Nf = 4 
case. The behavior is consistent with the hyperscaling Eq. The pseudoscalar mass 
dependence of the ratio at larger mass can be realized by the correction to the hyperscaling 
which may be different from one quantity to another. For /3 = 4 (open symbol), there is no 
fiat range without volume dependence. 
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FIG. 4. Dimensionless ratios F-j^/M-,^ and Mp/M-,^ as functions of oM,r for Nf = 12 at /3 = 3.7 
(filled symbol) and 4.0 (open symbol) for the two largest volumes. 

Similar observation can be made for the other ratio Mp/M^^ shown in the right panel 
of Fig. m Here, the flattening is observed for /3 = 3.7 again, but the range is wider than 
Ft^/Mt^. In this case, /3 = 4 shows the flattening, also. The difference of the constant is 
made possible due to a discretization effect. 

In the following, further detailed study using these data is performed. From the obser- 
vation here, we note that the region where hyperscaling is realized could be limited to the 
smaller masses. Further, the scaling range of the tt decay constant may be narrower at 
/3 = 3.7, and there may be no scaling range for the decay constant at /3 = 4. 

Existence of the scaling for at /3 = 3.7 and absence at /3 = 4 at the same aM^^ can 
be made possible if the M^^ in the physical unit is larger (thus the correction is no longer 
negligible) for /3 = 4, i.e. the lattice spacing decreases as /3 increases. In that case, the 
physical volume is smaller for /3 = 4, which gives a reason for the volume effect observed 
only for /3 = 4. 

To see if the inequality of the lattice spacing a(/3 = 3.7) > a(/3 = 4) holds, a matching 
between the two data sets of Mp/Af^- vs aAf^- is performed. We will not use F^/M^r for 
the matching due to the existence / nonexistence of the volume effect for /3 = 3.7 / 4. If 
an aMjr point which describes the same physics is found for each /3 separately, it can be 
regarded that the ratio of the lattice spacing as M^^ in physical units is the same for them. 
We exploit the mass dependence at the tail for this matching. To do that, flrst, the matching 
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FIG. 5. RMp/M-,^ is plotted against raM.^ after a crude matching with {R,r) = (1, 1) for /3 = 3.7 
and (1.035, 1.1) for /3 = 4. Legends are the same as Fig. [H 



of the overall factor absorbing the discretization error in the vertical direction needs to be 
performed, which can be done by introducing a factor R to multiply Mp/M^^ for /3 = 4 and to 
tune it to match the plateau to that of /3 = 3.7. Then, the remaining difference which may 
appear at the tail is absorbed by the horizontal matching factor r as r ■ M^^. Figure [5] is the 
same as Fig. HJ after a crude matching is done with {R,r) = (1.035, 1.1). In this particular 
definition of the lattice spacing and with a crude analysis, the ratio of the lattice spacing is 
a(/3 = 3.7)/(/3 = 4) ~ 1.1, and thus a{(3 = 3.7) > a{(3 = 4) holds. The fact that the lattice 
spacing decreases as (5 = G/g"^ increases is consistent with being in the asymptotically free 
domain, even if there is an IR fixed point in the beta function. 



IV. FINITE-SIZE HYPERSCALING 



A. Preliminary 



In the conformal window with finite masses and volume, the renormalization group anal- 
ysis tells us the scaling behavior for low-energy spectra which should obey the universal 
scaling relations S as 

Cp = LM, = f,{x), (10) 



For reviews, see, e.g., Refs. 
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FIG. 6. is plotted as a function of the scaling variable x for 7 = 0.1, 0.4, and 0.7 from left to 
right for Nf = 12 at /3 = 3.7. An alignment is seen for 7 ~ 0.4. 

where p distinguishes the bound state, p = tt or p in this study, or 

^F = LF^ = Mx). (11) 

The product of bound-state mass or decay constant and hnear system size falls into a function 
of a single scaling variable @ 

x = L-mJ+^, (12) 

where 7^, is the mass anomalous dimension at the IR fixed point. We call these scaling 
relations "finite-size hyperscaling" (FSHS). The forms of the scaling functions fp{x) are 
unknown in general. However, as the hyperscaling relation Eq. ([5]) must be reproduced in 
large volumes, the asymptotic form should be f{x) ~ x at large x. 

Now we examine whether our measurements of bound state mass and decay constant 
at different m/ and L obey FSHS. First, to visualize how the scaling works, we follow the 
analysis given, for example, in the model of the SU(3) with sextet fermions j47| and the 



SU(2) with adjoint fermions 48|. Panels in Fig. [6] show as functions oi x = L ■ mV^^'^'^^ 



for several values of 7. It is observed that the data points align at around 7 = 0.4, while 
they become scattered for 7 away from that value. This indicates a possible existence of 
FSHS with 7 ~ 0.4. A similar alignment is observed for in Fig. [7] as well. In this case, 
one finds the optimal scaling at around 7 ~ 0.5. 

These results are in contrast to our results for the Nj = 4 system with the same lattice 
action with /3 = 3.7, in which the chiral symmetry is spontaneously broken. It is found 
in Fig. |8]that the alignment is observed at 7 = 1, which is interpreted as a realization of 



^ We adopt this definition, which corresponds to x^^'^* of Ref. 



451 ] ■ to make its asymptotic form simpler. 
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FIG. 7. as a function of x for 7 = 0.2, 0.5, 0.8 from left to right for Nf = 12 at f3 = 3.7. An 
alignment is seen for 7 ~ 0.5. 
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FIG. 8. is plotted as a function of the scaling variable x for 7 = 0.5, 1.0, and 1.5 from left to 
right for Nf = 4 at /3 = 3.7, where the S^SB occurs. An alignment found at 7 = 1 is consistent 
with Eq. 



Eq. ([6]). The vr decay constant does not exhibit aUgnment at any value of 7 allowed for the 

(Fig. ED- 



unitarity requirement < 7 < 2 
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B. Quantitative analysis 



To quantify the "alignment" we introduce an evaluation function P(7) for an observable 
p as follows. Suppose is a data point of the measured observable p at 
and 5^^ is the error of j labels distinction of parameters L and m/. Let be a subset 
of data points {{xk^^k)} from which we construct a function f^^\x) which represents the 
subset of data. Then, the evaluation function is defined as 

^(7) = T7>.>. ^r-^ . (13) 
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FIG. 9. as a function of x for 7 = 0, 1, and 2 from left to right for A'^^ = 4 at /3 = 3.7. No 
alignment is seen in this case. 

where L runs through all the lattice sizes we have, the sum over j is taken for a set of data 
points which do not belong to Kl which includes all the data obtained on the lattice with 
size L. J\f denotes the total number of summation. Here, we choose for the function f^^^^ a 
linear interpolation of the data points of the fixed lattice size L for simplicity, which should 
be a good approximation of ^ for large x. 

This evaluation function takes a smaller value when the data points are more closely 
collapsed to the line f^^^^ and thus provides a measure of the alignment. P(7) varies as the 
choice of parameter 7 and should show a minimum at a certain value of 7 when the optimal 
alignment of data is achieved. We take it as the optimal value of 7. 

We then estimate the uncertainty of the optimal 7 by properly taking account of the 
statistical fluctuation of as well as its effect to the line f^^^\ For this purpose, we employ 
the parametric bootstrap method, in which the data point is simulated by a random sample 
generated by Gaussian distribution with the mean ^-^ and the standard deviation 6$,^ . The 
distribution of 7 is thus obtained for a large number of these samples, from which the 
variance of 7 is estimated. The systematic error associated with the interpolation will be 
estimated by choosing different functional form with linear or quadratic splines as will be 
discussed subsequently. 



Our method for lattice study is similar to those used in the literature |32|, |47|, |52| , based 



on the study of the critical phenomena in a finite-volume system [53| l|. We incorporated the 
uncertainty of data points as the weights in the evaluation function so that it is normalized 
to one when the distance between the data points and the interpolated line is equal to the 
^ After submitting this paper, we were informed that a similar method had been considered in the Utera- 
ture 5J, |55| . We thank David Schaich for the information. 
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FIG. 10. (Left): The 7 dependence of the evaluation function P for M^- at /3 = 3.7 is shown. The 
vertical axis shows the central values of P as a function of 7. Each curve shows the results of -P(7) 
with the corresponding range. (Right): The results of P{^) using data sets with two different 
volumes as L = 18, 24, L = 18, 30, and L = 24, 30, respectively, in which the whole x range is 
considered. 



standard deviation of the data point. The unweighted version of the evaluation function has 
also been examined in our analysis, which resulted in values consistent with the optimal 7. 

In the evaluation function Eq. ( IT3l) . the data points need to be taken for a range of 
X = L ■ rri^j.^^^^"'^ in which there is an overlap of available data for all volumes, L = 18, 24, 
and 30. The maximum value oi mf is chosen so that ~ 1 is satisfied, and the minimum 
is chosen so that the finite-volume effect on the mass of the bound state is not too large. 
Therefore, for all the values of 7 to be examined, we consider the range of as follows: 
The minimum Xmin should be such that frif = 0.04 for the largest volume L = 30 at /3 = 3.7, 
OT rrif = 0.05 at /3 = 4, and the maximum Xmax should be such that mj = 0.2 for the 
smallest volume L = 18 at /3 = 3.7, or mj = 0.24 at /3 = 4. Around the optimal value of 
7, we have 12 data points for /3 = 3.7 and 11 for /3 = 4 within the range [xmin? a^max]- Note, 
however, we may need to incorporate some neighboring data outside this range to obtain 
the interpolated value f^^\x) by the spline functions. 

The line marked by "all" in Fig. [10] plots the values of the evaluation function using all 
data in Xmin < a; < Xmax for /3 = 3.7. A clear minimum is observed, at which the optimal 
alignment of the data is achieved. We repeat this analysis for each observable and at each 
value of (3. The results are tabulated in the column labeled by "all" of Table I VI II Figures 
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FIG. 11. The data for M-j^ (left), F^^ (center), and Mp (right) are plotted as a function of x at 
/3 = 3.7. The vertical axis shows the values for ^. The horizontal axis shows the x = L - rnj^^~^^^ 
where the values of 7 are 0.436, 0.513, and 0.476 for Ad-j^^ F-^-, and Mp, respectively. 
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FIG. 12. The data for (left), (center), and Mp (right) are plotted as a function of x at 
/3 = 4. The vertical axis shows the value for ^. The horizontal axis shows the x = L ■ rri'j^^^"'^ 
where the values of 7 are 0.408, 0.560, and 0.479 for M^, F^^, and Mp, respectively. 



[TT] and [12] plot the values of C, for each case as a function of x with the optimal 7. 

Let us now recall that from the naive analysis of the ratio in Sec. Illlt the scaling region 
may be restricted to a range of smaller masses. To capture such an effect, we need to 
study systematically how the range of fermion mass affects the evaluation function and the 
resulting optimal 7. In this regard, we define a window of the parameters x and L to which 
the data are restricted in evaluating the evaluation function and see if the results change 
against the choice of the window. 

We consider a window in the x direction which has a span Ax/2, where Ax = x^g^^ — x^i^ 
and slide it with an interval Ax/4. There are three choices of windows denoted by range 1, 
2, and 3, which are [xmm, a:min + Ax/2], [xmin + Ax/4, x^ax- Ax/4], and [xmax - Ax/2, x^ax], 
respectively. In these windows, the data are taken for all three volumes. For the L direction 
we have three windows: The first window consists of data for L = 18 and 24, the second 
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TABLE VII. Summary of the optimal values of 7. See the text for details. 



quantity 


/9 


all 




X 






L 




range 1 


range 2 


range 3 


(18,24) 


(18,30) 


(24,30) 




3.7 


0.434(4) 


0.425(9) 


0.436(6) 


0.437(4) 


0.438(6) 


0.433(4) 


0.429(8) 




4 


0.414(5) 


0.420(7) 


0.418(6) 


0.411(5) 


0.397(7) 


0.414(4) 


0.447(9) 




3.7 


0.516(12) 


0.481(19) 


0.512(19) 


0.544(14) 


0.526(18) 


0.514(11) 


0.505(24) 




4 


0.580(15) 


0.552(21) 


0.602(20) 


0.605(19) 


0.544(27) 


0.577(14) 


0.645(32) 


Mp 


3.7 


0.459(8) 


0.411(17) 


0.461(10) 


0.473(8) 


0.491(15) 


0.457(8) 


0.414(18) 


Mp 


4 


0.460(9) 


0.458(13) 


0.455(14) 


0.460(8) 


0.457(16) 


0.459(8) 


0.463(15) 



window for L = 18 and 30, and the third window for L = 24 and 30. In these windows in L 
directions, the whole x range [xmin,a:max] is considered. 

The evaluation function P(7) of the vr mass for each window is also plotted in Fig. [TOl 
It is noted that the value of -P(7) is 0{1) at the minimum. The optimal values of 7 which 
minimize -P(7) are summarized in Table IViTt where the numerals in the parentheses include 
the statistical and systematic errors. The systematic error due to the ambiguity of the 
interpolation is estimated by the difference of the optimal 7's obtained with linear and 
quadratic spline interpolations. The comparison of these -P(7) is shown in Fig. [131 The 
minima for the quadratic spline interpolation appear approximately at the same place as 
those for the linear one. The systematic error is thus small, and it is found to be always 
smaller than the statistical error in the analysis with the whole data points. 

Let us first look at the case of M^r at /3 = 3.7 in Table IVTll We do not observe the window 
dependence beyond the error bars, and thus is well-described by a function of a single 
variable x. If there is IR conformality, the nonuniversal correction to the hyperscaling is 
negligible at this precision. Then, from the fact that in Sec. IIIIl the scaling is observed in 
the small mass range for the ratio Mp/M^ at /3 = 3.7, it is suggested that there should be 
certain window dependence of 7 from Mp. As 7(iVf7r) = 0.434(4) and 7(-F7r) = 0.459(8) at 
/3 = 3.7, if one restricts the mass range for Mp to the smaller side, then the value of 7(Mp) 
should get closer to that of 7(71/7^). This is actually the case, as is observed from Table lylTl 
in which 7(Mp) reduces for smaller mass range (toward range 1) and larger volume (toward 
L = 24, 30). 
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FIG. 13. The 7 dependence of the evaluation function P for Mj^, Fj^, and Mp at 13 = 3.7 is 
plotted. The vertical axis shows the values of P at each of 7 where the three volumes and full 
range of x for the data are considered. The solid and dashed curves show the results of P{'j) with 
the interpolation functions f{x) by the linear and quadratic functions, respectively. 



However, such trends are not observed for the Mp at /3 = 4, where one expects the similar 
X and L range dependence. As the number of samples have gotten reduced for the fixed 
range analysis, a statistical instability might have spoiled the result. 

The similar trend for the x-range dependence as for Mp at /3 = 3.7 is observed for at 
/3 = 3.7, too. The direction of the movement is correct, but it does not get close enough to 
the value of 7(M.n-). Moreover, the L range dependence is too weak to conclude that it will get 
close to '-/{Mj^). These results may be understood from the fact that in Sec. IIIH the scaling 
is observed only in the very small mass range. For at /3 = 4, the L dependence appears to 
be opposite to the expectation, which can be understood as the result of unobserved scaling 
in the analysis in Sec. IIIII 

As we cannot completely resolve these trends in the mass dependence, we regard these 
variations of 7 with respect to the change of the window as the systematic error on the 
central value of 7 obtained with "all" data. We put the asymmetric error for both x and 
L directions separately estimated by the maximum variations from the central value. The 
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FIG. 14. The results of the values of 7 for three observables at two /3 are summarized, where the 
statistical and systematic errors are added in quadrature. All the results are consistent with each 
other within 2a level, except for the 7 from at /3 = 4, where the scaling region is suspected to 
be outside of the parameter range we have. See the text. 



results read 



7 



0.434(4) (l^)(l?o) 
0.516(12)(l}?)(ti) 

0.459(8) (11) (1^^) 
0.414(6) (1?^)(+D 

29 J 



0.580(15)(1^?)( 
0.460(8)(1^)( 



for 
for 

for Mp 

for 

for 

for M„ 



> for /3 = 3.7, 



(14) 



> for /3 



where the numerals in parentheses are the combined errors of the statistics and the in- 
terpolation, the systematic uncertainties due to the dependence on the volume, and the 
systematic errors due to the dependence on the x range, respectively. The results with all 
the errors added in quadrature are summarized in Fig. [TH All the results are consistent 
with each other within 2a level, except for 7 from at /3 = 4 for which the scaling region 
was suspected to be outside of the parameter range we have examined. (See Sec. IIIII ) 

From these analyses, we conclude that our data for the Nf = 12 theory by the lattice 
simulations are reasonably consistent with the FSHS and, thus, with the IR conformality, 
if we exclude Fj^ at P = 4 from the analysis. The resulting mass anomalous dimensions 
from different quantities and two lattice spacings are also reasonably consistent. We quote 
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0.4 ^ 7* ~ 0-5 for the value of the mass anomafous dimension at the IR fixed point. 



V. NONUNIVERSAL CORRECTIONS TO FINITE-SIZE HYPERSCALING 

We found in the previous sections that our data is consistent with the FSHS once the 
systematic effect due to the hmited parameter range is taken into account. The analyses 
also suggested that the nonuniversal correction to the hyperscaling could be important. In 
this section, we try to test two plausible models for the correction. To do this test, we 
need to fix the term for the universal scaling in the following. Therefore, the approach loses 
generality that the analysis in the previous section had. Thus, the result here is not going 
to be the main result in this paper, but it still provides useful information. 

For this purpose, we use the following formulae for the fit: 

e = Co + ciLmy(^+^) ■■■fit A, (15) 
e = Co + ciLm)^'-^^^^ + csLm^ ■ ■ ■ fit B. (16) 

The fit A uses a naive fit form based on the hyperscaling relation which is described by the 
function form of f{x) = Co + CiX with x = Lni'^^ This formula is motivated from the 
results obtained in Fig. [Ill since the clear linearity of the data for large x can be found near 
the optimal value of 7. The second formula for fit B includes a mass correction term. 

As discussed in the previous sections, there may exist some corrections beyond the hy- 
perscaling relations in the region we simulated, so we try to include such contributions. In 



general, finite- volume corrections exist, and the form is given a la Fisher [56|, l57|]. We do not 
take into account these effects, however, as the analysis in the previous sections indicates 
the mass correction is more important for the parameter space we simulated. 

Among various possible choices, we take a = (3 — 27)/(l + 7), which is inspired by the 
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. We also consider 



analytic expression of the solution of the SD equation given in Ref. 
a = 2 case, which is regarded as the small mass correction caused by the explicit chiral- 
symmetry-breaking effects or due to the lattice discretization artifact. It is noted that, in 
both cases, the fit function cannot be described by a single scaling variable x = Lm^j^^^'^\ 
We denote these fit functions with a = (3 — 27)/(l + 7) and a = 2 as fit Bl and fit B2, 
respectively. All the details of the fit results are shown in Appendix [B], which also includes 



other fits, for example, the case of a 



3l| . We only give a digest of the results. 
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TABLE VIII. The fit results of finite-size conformal hypothesis at /3 = 3.7(left) and /? = 4(right). 
The values in the brackets mean the input in the fit. 





7 


a 


x7dof 


dof 




7 


a 


xVdof 


dof 


fit A 
fit Bl 
fit B2 


0.449(3) 
0.411(9) 
0.423(7) 


(3-27) 
(1+7) 

[2] 


4.52 
1.23 
1.15 


47 
44 
44 


fit A 
fit Bl 
fit B2 


0.430(3) 
0.461(18) 
0.453(11) 


(3-27) 
(1+7) 

[2] 


6.78 
1.86 
2.00 


44 
41 
41 



We perform simultaneous fits using M^., Mp and F^, with common 7 for each fit ansatz: 
fit A, Bl, B2, and for each /3 separately. In these fit analyses, we assume for simplicity 
that possible correlation between the data of ,^7^, C,f and is absent. Note that, here, 
we include the data at /3 = 4, which were excluded in the analysis in the previous 
section due to the absence of manifest hyperscaling (Fig. H]). This is because the mass 
correction in the hyperscaling analysis, which is included here but not there, could amend 
the possible disturbance to the hyperscaling from the mass correction. To avoid the large 
finite-size effects, the data used in the fits are restricted to those which satisfy = LM^^ > 9 
determined by trial and error. Figures [15] and [16] show the fit B2 as examples. Table IVIIII 
shows the resulting fitting parameters. 

The x^/dof is large ^ 5 for the fit A, while it is reasonable (x^/dof < 2) for the fit Bl and 
B2. This implies that the correction is necessary, and both types of mass correction work. 
Let us look into the fit B2 at /3 = 3.7, for example. The ratio of the terms proportional to 
Ci (leading scaling) and C2 (mass correction) at the heaviest point {nif = 0.2) is —0.008(4) 
(M^), -0.034(6) (Mp) or -0.067(7) (F^). The small correction for is consistent with the 
fact that the value of 7(^,1) was stable against the x and L range used in the analysis in 
Sec. IIVI Furthermore, the order of the size of the correction for M^, Mp, and F^- is consistent 
with what we observed in Sees. IIIII and IIV[ 

We obtain consistent 7 for all the fit Bl and B2. The resulting values are in the range 
7 ~ 0.41 — 0.46, which is consistent with the main result presented in Sec.[IVl 0.4 ^ 7^, ^ 0.5. 
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FIG. 15. Spectra ^^(left), ^^^(center), and ,^p(right) as a function of m/ at /3 = 3.7. For simplicity, 
we only show two fit results of fit A and fit B2, by the solid and dotted curves, respectively. The 
data with empty symbols are not used in the fit. 




FIG. 16. Spectra ^^^(left), .^i? (center), and i^p(right) as a function of at /3 = 4. For simplicity, 
we only show two fit results of fit A and fit B2, by the solid and dotted curves, respectively. The 
data with empty symbols are not used in the fit. 

VI. CHIRAL PERTURBATION THEORY ANALYSIS 

A counter scenario to the conformal hyperscaling is the S^SB. Here, we give a digest of 
the analysis, the detail of which is described in Appendix O 

As the current data set of /3 = 3.7 simulation has physically lighter quark mass, the 
property of the chiral limit is captured better than those at /3 = 4. Therefore, we focus on 
the former. First, we take the infinite volume limit of and F^^ using a formula inspired by 
ChPT or finite-size effect of the multiparticle state a la Liischer. Then, the vr decay constant 
and mass squared are fit with polynomial functions in m^. Logarithmic mass dependence of 
these quantities would emerge at the chiral regime and play an important role there. As it 
turns out, however, our mass range is far from it. 

Several fitting ranges are examined for with second order polynomial in m/. With 
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fixed minimum m™™ = 0.04, the maximum is changed as m™*^^ = 0.12, 0.1, and 0.08. 
Only the m™*^^ = 0.08 fit is consistent with the vanishing intercept, while others end up 
with having negative intercept. A fit with the intercept constrained to be zero gives good 
X^/dof = 0.88 for m'j^^ = 0.08, though including heavier mass gives unacceptably large 
xVdof > 8. 

Similar exercises are done for the F^, where the intercept is always positive, and for 
j^max = Q I g^Yid 0.08, reasonable x^/dof yields x^/dof = 0.37 and 0.29, respectively. 

The 71 mass squared and decay constant behave in the way reasonably consistent with the 
second-order polynomial, when the mass range is restricted to smaller side 0.04 < nif < 0.08 
where we have four mass points. The tt mass vanishes in the chiral limit, and the decay 
constant is nonzero. However, we are not able to conclude these results are consistent with 
SxSB. The problem can be appreciated if the expansion parameter in ChPT is evaluated. 
Jsing the value of in the chiral limit F^^ = 0.0190(52), the natural expansion parameter 
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X = Nf{ ) ^39, (17) 



- |l5| of ChPT in the large Nf theory as given in Eq. ([T]) is 

2 

47iFjV2j 

which has been evaluated at the lightest point rnf = 0.04 on L = 30 lattice, where = 
0.3028(14). For the expansion to be consistent, it should satisfy X < 1. With this large 
expansion parameter at the smallest mass in the deposit, we are not in the position to judge 
if the decay constant is really nonzero at the chiral limit from the extrapolation performed 
here. 



VII. SUMMARY AND DISCUSSION 

We have studied the SU(3) gauge theory with 12-fiavor fermions in the fundamental 
representation of the gauge group. We tried to determine if the massless theory undergoes 
SxSB or exhibits IR conformality due to the existence of an IR fixed point of the beta 
function, through the fermion mass rrif dependence of the spectrum. A type of HISQ action 
was adopted to maximally suppress the staggered flavor-symmetry violation as well as other 
discretization errors. Three volumes with fixed aspect ratio and various rrif were examined. 
Simulations were repeated for two values of the bare gauge coupling, corresponding to two 
different lattice spacings. We observed that the lattice spacing became finer when the bare 
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gauge coupling was made small, which is expected when the simulation points are in the 
asymptotically free domain. 

A primary analysis of the masses of pseudoscalar and vector channel and the pseudoscalar 
decay constant revealed a range of small rrif where the ratios of composite masses and 
decay constant, Ft^/Mt^ and Mp/M.,^, are independent of rrif. This is consistent with the 
hyperscaling characteristic to the IR conformality and was observed for all the quantities 
except for the one which involves the decay constant at a finer lattice. 

Further detailed study adopting a finite-size scaling have shown that our data are reason- 
ably consistent with the FSHS, where the product of linear system size and the composite 
masses (or decay constant) fall into a function of a universal scaling variable composed of 
the rrif, linear system size L, and mass anomalous dimension 7=,, at the IR fixed point at 
low energy. The resulting 7^,'s obtained from the introduced evaluation function, -P(7), of 
the scaling were reasonably consistent with each other for three observables and two lattice 
spacings when only the decay constant at finer lattice was excluded as before. We conclude 
the mass anomalous dimension is in the range 0.4 ^ 7^, ^0.5 at IR fixed point if it exists. 

Existence of the nonuniversal correction indicated by the FSHS motivated a study of 
global fit with models assuming some corrections. By adding a correction term to the 
lowest-order universal term, the global fit works well with a reasonable x^? though possible 
correlation was assumed to be absent. Among the tested, the one with the power of the 
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and another with 0((ma)^) discretization error 



fermion mass fixed a la ladder SD study 
resulted in 7* consistent with FSHS. 

In the test of our data against the S^SB scenario, the ruf dependence of the vr mass 
squared and decay constant appeared to be consistent with the tree-level mass dependence 
to second order from chiral perturbation theory, if the data are restricted to sufficiently small 
rrif. The decay constant remains nonzero at the chiral limit with this analysis. It turns out, 
however, the natural expansion parameter of ChPT estimated using the value of the decay 
constant at the chiral limit extracted in this analysis is very large X even at the smallest tt 
mass simulated. Therefore, we cannot judge if the decay constant is really nonzero at the 
chiral limit from the extrapolation performed in this study. 

A possibility of S^SB in A'^^ = 12 is not excluded yet. Further efforts would be required 
to arrive at a decisive conclusion. One possible direction is to use the brute force (larger 
volume and lighter mass) calculations, which seems very hard. Another may be to gather 
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more information about the spectrum, for example, the mass of glueball and flavor singlet 
composite state, which could have distinctive signature across the phase boundary. 

However, we have confirmed that the general property of the spectrum in Nf = 12 theory 
as functions of is completely different from the Nf = 4 theory [s^], where we found that 
the ratio Fj^/Mt^ gets divergent towards the chiral limit in accord with the ordinary QCD but 
in sharp contrast to the conformal behavior we have observed for Nj = 12. Although the 
analyses support the conformal scenario, even if the Nf = 12 theory breaks chiral symmetry, 
the breaking scale is very small, so it must be very close to the boundary of the conformal 
phase transition point. In this case, our result of 7m obtained with assuming the IR fixed 
point can be regarded as an approximate mass anomalous dimension at the walking regime. 

There have already been several studies focusing on the conformality of the Nf = 12 
SU(3) gauge theory by investigating the masses of composite states of the mass-deformed 
theory: 

Ref. 2J] performed a large-scale simulation using the stout-smeared staggered fermions 
at a single lattice spacing. They performed the global fits of ChPT and conformal scenario. 
ChPT fit (without logarithmic terms as in ours), which resulted in a nonzero pion decay 
constant and chiral condensate, was favored over conformal scenario by examining the x^/ dof 
of the fits. However, the expansion parameter at their lightest quark mass reads X ~ 34, 
which is as large as ours [Eq. ( fT7j) ]. We note that we have concluded that for such a large 
X, ChPT is not self-consistent, and hence it is difficult to conclude what is really happening 
at the chiral limit. 

As to their conformal fit, it is equivalent to our fit A in Eq. ( fT5|) . by which we obtained 
a x^/dof similar to theirs. (Their estimate of the effective mass anomalous dimension reads 
7* = 0.395(52), which is consistent with ours.) Although it is not our main analysis, we 
introduced a mass correction to the conformal type fit (fit Bl or B2), which was required 
by our data, and obtained a reasonable x^/dof. Looking into more details, we observe a 
difference in the quark- mass dependence of pion decay constant. Their decay constant is 



perfectly fitted with the linear function in m/ (Fig. 1 middle in Ref. 2J]), while ours exhibits 



a curvature (Fig. [15] middle or Fig. [2T!) . The origin of different behavior is so far unknown. 



There have also been studies analyzing the data of Ref. 2J] for tests of the conformal 



hypothesis. Ref. 3l| performed fits with the single power corresponding to the fit A in 



Eq. f|T5|) without the constant term and examined the dependence of as function of input 
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7- The optimal vale of 7 which minimizes reads 7 ~ 0.4 for vr and p masses, which 
is consistent with our result. However, they obtained 7 ~ 0.25 for pion decay constant, 
while our ^{F^^) obtained without fixing the function form in Sec. [IV] tends to be even larger 
than 7(M7r). This inconsistency likely originates from the aforementioned difference in mass 
dependence. They tried the global fit to various observables with a common 7 introducing 
a volume correction term to every observable and a mass correction term to the pion decay 
constant. By doing that, a model dependence is introduced, but our result is consistent with 
their result: 7 = 0.403(13) with a reasonable x^/dof = 42/44. 



Ref. 
Ref. 



32| tried a similar analysis to what is presented in Sec. lIVBt using the data of 



24j . and obtained 7 ~ 0.35. The error is large (±0.23 for pion mass), thus consistent 



with our result. 
Ref. 



21| performed a pure power-law fit to vr and p masses on a single volume but over 
two lattice spacings. The results are 7 ~ 0.52 — 0.64, which is only a bit higher than what 
we obtained. 

To conclude, the (effective) 7* did not appear as large as the walking theory should exhibit. 
By this, the hunt for the realistic walking theory should anyway alter the direction of the 
number of flavors towards smaller. The Nf = 8 theory, which could have large anomalous 
dimension 7* ~ 1 as a candidate of the walking technicolor model, is under investigation 



with exactly the same setting with HISQ action 
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Appendix A: Staggered flavor symmetry 



We show some other staggered meson spectra e.g. non-NG pseudoscalar (7574 ® ^5^4) 
denoted as SC and a vector channel (74 (8 ^i) denoted as VT in comparison with the cor- 
responding (NG) 71, PS: (75 ® ^5) and the vector, PV: (7i74 ® ^4^4), which are used in the 
analysis given in the main text. This is to see a staggered flavor-symmetry-breaking effect 
in our HISQ simulation. The results are shown in Fig. [171 As we expected, the two masses 
for each pseudoscalar and vector meson channel are almost degenerate. 





FIG. 17. Result of mass spectra of four operators at /? = 3.7 (left panel), 4 (right panel) on 
(L,r) = (30,40). 



Appendix B: Details of fit analysis for conformal hypothesis with corrections 

Numerical detail of fit results on the conformal hypothesis are given here. In the conformal 
hypothesis with a finite volume, we make an attempt to use the fit functions given in 
Eq. [151 We fix the value of the exponent a to a certain value or to relate with 7 in the fit 
since having two exponents free makes the fit unstable. We consider three possible cases: 
a = (3 — 27)/(l -|- 7), 2, and 1; we denote these fit functions as fit Bl, fit B2, and fit 
B3, respectively. We carry out simultaneous fits using all the data for M^^, F^^, and Mp 
with common anomalous dimension 7 and/or the parameter a where possible correlation is 
assumed to be absent. We use same data set for the fit as given in Sec. [Vl The fit results are 
shown in Table [IXl and [Xl for /3 = 3.7 and /3 = 4. As already discussed in Sec. [Vl additional 
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TABLE IX. The fit results of finite-size conformal fiypothesis at /3 = 3.7. 



fit A 


Co 


Cl 


fit Bl 


Co 


Cl 


C2 




-0.029(28) 


2.793(5) 




0.095(30) 


2.924(39) 


-0.227(83) 




0.147(12) 


0.544(3) 




0.627(10) 


-0.230(27) 


0.037(17) 




0.363(58) 


3.325(10) 


if 


0.146(77) 


3.634(52) 


-0.74(14) 



7 = 0.449(3), xVdof=4.52, dof=47 7 = 0.411(9), xVdof=1.23, dof=44 



fit B2 Co Cl C2 fit B3 cq ci C2 





0.089(30) 


2.860(26) 


-0.181(96) 




0.095(30) 


3.53(25) 


-0.81(27) 




0.045(16) 


0.600(7) 


-0.325(33) 


Cf 


0.030(17) 


0.829(59) 


-0.356(67) 


ep 


0.147(75) 


3.524(37) 


-0.97(16) 




0.149(79) 


4.56(32) 


-1.45(36) 




7 = 0.423(7), 


xVdof=1.15, 


dof=44 


7 


= 0.356(22), 


xVdof=1.45, 


dof=44 



correction terms improve x^/dof for both cases of /3 = 3.7 and /3 = 4. In both cases of fit 
Bl and fit B2, the values of 7 are consistent with the value from the analysis of the FSHS 
test given in Sec. llVi On the other hand, the result of fit B3 gives a slightly smaller value 
of the 7, but with larger errors. 



Appendix C: Details of chiral perturbation theory analysis 

In order to give a fair comparison whether or not the chiral symmetry is spontaneously 
broken, we carry out the fit based on the ChPT hypothesis by using our data. 



1. Finite-size dependence of physical quantities 
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Here, we study the finite volume effects for M^^ and F^. To obtain t 
Ft,- in the infinite volume limit, we use the following ChPT-inspired 

M^{L) -M^ = CM^ 



le values of and 



F^{L)-F^ = cf^ 



60| fit functions: 

(Cl) 
(C2) 



31 



TABLE X. The fit results of finite-size conformal hypothesis at /3 = 4. 



fit A 

ill Jr\. 


Co 




fit R1 


Co 


Cl 


C2 




-0.202(31) 


2.662(5) 




-0.083(33) 


2.486(74) 


0.33(12) 




0.212(15) 


0.464(3) 




-0.040(25) 


0.531(15) 


0.166(34) 




0.228(60) 


3.036(9) 




-0.146(84) 


3.013(85) 


-0.03(16) 


7 = 


0.430(3), xVdof=6.78, 


dof=44 


7 


= 0.461(18), 


xVdof=1.86, 


dof=41 



fit B2 Co Cl C2 fit B3 cq ci C2 





-0.096(33) 


2.551(34) 


0.45(12) 




-0.074(33) 2.30(21) 


0.44(24) 


^F 


-0.026(25) 


0.518(9) 


-0.291(36) 




-0.044(25) 0.564(54) 


-0.147(62) 




-0.123(81) 


3.039(42) 


-0.16(16) 




-0.153(86) 2.92(27) 


0.09(31) 


7 


= 0.453(11), 


xVdof=2.00, 


dof=41 


7 


= 0.476(39), xVdof=1.88, 


dof=41 



where cm^ and cp^ are the fit parameters. We carry out the simultaneous fit for each fermion 
mass by using three data points on L = 18,24, and 30 at /3 = 3.7 and /3 = 4, where we 
assume for simpHcity that possible correlation between the data of M^. and is absent. 
The fit results are shown in Figs. [18] and [191 As a result, in the entire fermion mass region, 
our data is well-fitted, and the value of x^/dof for each parameter is 0{1). Also, one can 
find that the difference between the value of L = 30 data and that in the infinite volume 
limit is negligible. Thus, we use the data on L = 30 to analyze the chiral behaviors of both 
and F^^ hereafter. 

2. ChPT fit analysis 

We analyze the fit using the second-order polynomial function as 

h{'m f) = Cq + Cim f + C2m^f (C3) 

where Co,i,2 are fit parameters. Using this simple polynomial function, we carry out the 
fits for and by the function h{mf), varying the fit range of the fermion mass from 
rrif = 0.04 to rrif = 0.12. We denote the fit range of fermion mass as [m™™,m™'^^]. As the 
current data set of /3 = 3.7 simulation has physically lighter quark mass, the property of the 
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FIG. 18. The results of the finite- volume scaling fit for and i^^ at f3 = 3.7 using the infinite 
volume extrapolation formula in Eqs. (jCip and ()C2p . 
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FIG. 19. The results of the finite volume scaling fit for M^- and F,r at /3 = 4 using the infinite 
volume extrapolation formula in Eqs. (jCip and ()C2p . 

chiral limit is captured better for /3 = 3.7 data than for /3 = 4 data. Therefore, we focus on 
the former here. 

The fit results for are shown in Fig. |20] and Table IXII Among these fits, the best-fit 
values of intercept cq is consistent with zero for the fit with mj-"-^ = 0.08, while we obtained 
negative values of Cq for others. In Table \XI\ we also show the fit results with cq being fixed 
to zero. As we expect, the values of x^/dof become much larger than unity except for the 
case with m™^^ = 0.08. Furthermore, the contribution of the higher-order term with C2 is 
not small enough even in the fit result using the data with the smallest mass range. It may 
indicate that the fermion mass in the data we have is too heavy to take a reliable chiral 
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extrapolation with ChPT-like formula. 

In Fig. [21] and Table |XII| we show the fit results for the case of F-,^. The same fit function 
and fit ranges as in the case of are used for the fits. Nonzero value of in the chiral 
limit is an important signal of the S^SB. Results show 0{1) value of x^/dof for the fit with 
smaller mass range, and we obtain tiny but nonzero value of in the chiral limit. 

If we adopt the fit result with the mass range [0.04, 0.08] at /3 = 3.7, which is expected to 
be the most reliable result among those fits, it might look like they are consistent with the 
hypothesis of S^SB scenario because the result shows that the pseudoscalar mass squared is 
going to zero along the tree-level ChPT, while the decay constant remain nonzero (though 
it is tiny) in the chiral limit. However, to be able to conclude that it is consistent with 
the ChPT fit, one has to check that the expansion parameter of the perturbation is in the 
legitimate region. Using the value of in the chiral limit, F^^ = 0.0190(52), the value of the 
expansion parameter is estimated as A" ^ 39. Here, we evaluated it using the data with the 
lightest mass, rrif = 0.04, on L = 30 lattice, where M^, = 0.3028(14). Since the expansion 
parameter is much larger than one even in the lightest mass region of the fit range, we are 
not in the position to judge if the decay constant is really nonzero in the chiral limit from 
the extrapolation performed here. Therefore, investigation with larger volume and smaller 
masses is needed to draw more definite conclusions regarding the test of the ChPT scenario. 
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FIG. 20. The several fit results for at /? = 3.7 using the data on L = 30. 
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TABLE XL The fit results on at (3 


= 3.7 using 


the data on L = 30. 


The values in 


the bracket 


mean the inputs 


in the fit. 










fit range 


Co 


Cl 




xVdof 


dof 


[0.04, 0.08] 


-0.0057(91) 


1.82(32) 


15.2(2.6) 


1.35 


1 




[0] 


1.62(3) 


16.76(45) 


0.88 


2 


[0.04,0.1] 


-0.0209(48) 


2.37(15) 


10.6(1.1) 


2.59 


2 




[0] 


1.729(21) 


14.99(25) 


8.33 


3 


[0.04,0.12] 


-0.0183(31) 


2.28(87) 


11.21(55) 


1.90 


3 




[0] 


1.780(17) 


14.28(17) 


10.29 


4 



TABLE XII. The fit results on i<V at /? = 3.7 using the data on L = 30. 



fit range 


Co 


Cl 


C2 


xVdof 


dof 


[0.04, 0.08] 


0.0190(52) 


1.21(18) 


-2.2(1.5) 


0.29 


1 


[0.04,0.1] 


0.0162(30) 


1.31(85) 


-3.01(58) 


0.37 


2 


[0.04,0.12] 


0.0231(18) 


1.093(48) 


-1.51(29) 


3.30 


3 
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